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C3 • 1 INTRODUCTION 



Observational evidences accumulated over the past few years in- 
dicate that the physical properties of supermassive black holes 
[SBHs, hereafter] residin g at the centr e of most, if not all, 
spheroidal galaxies (Richstone et al. 1998), correlate with those of 
the host galaxy iMagorrian et alJll998tlFerrarese & MerritdEoOOt 
lGebhardtetail2 000) and, possibly, of the host dark matter (DM) 
halo lFemu~ese 12002). Although it is not clear whether the ob- 
served correlation between halo and BH masses is genuine or sim- 
ply reflects the fact that massive haloes preferentially host massive 
spheroids (Wvithe & Loeb 2005a), the strong connection between 
the BH properties and the gravitational potential wells that host 
them suggests a link between the assembly of BHs and the evolu- 
tion of galaxy spheroids. This co-evolution is expected within the 
frame of hierarchical models of structure formation, like the pop- 
ular cold dark matter (CDM) cosmogonies. Since it is generally 
accepted that the quasar activity is powered by SBHs, several mod- 
els have been proposed to trace the assembly of SBHs during the 
hierarchical build-up of their host haloes using numerical, analytic 
and semi-analytic methods. In this work we focus on the latter two 
approaches that consist of parameterizing in term of simple analytic 
models the complex physics of the galaxy formation process, while 
the evolution of DM haloes can be either followed using numerical 



ABSTRACT 

We investigate the ability of hierarchical models of QSO formation and evolution to 
match the observed luminosity, number counts and spatial clustering of quasars at redshift 
Z < 2. These models assume that the QSO emission is triggered by galaxy mergers, that the 
mass of the central black hole correlates with halo properties and that quasars shine at their 
Eddington luminosity except, perhaps, during the very early stages of evolution. We find that 
models based on simple analytic approximations successfully reproduce the observed fi-band 
QSO luminosity function at all redshifts, provided that some mechanisms is advocated to 
quench mass accretion within haloes larger than ~ 1O 13 M that host bright quasars. These 
models also match the observed strength of QSO clustering at z ~ 0.8. At larger redshifts, 
however, they underpredict the QSO biasing which, instead, is correctly reproduced by semi- 
analytic models in which the halo merger history and associated BHs are followed by Monte 
Carlo realizations of the merger hierarchy. We show that the disagreement between the lumi- 
nosity function predicted by semi-analytic models and observations can be ascribed to the use 
of B-band data, which are a biased tracer of the quasar population, due to obscuration. 
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techniques, in the so called semi-analytic approach, or simplified 
analytic prescriptions, which constitute a full analytic model. 

Over the years many different analytic (see, e.g., 
Efstathiou & Reesfll988t lHaehnelt & Reeslll993l: lHaiman & Loebl 
19981 IPercival & Milled Il999t lHaiman & Menod |200(j 
Martini & Weinberg 2001; Hatziminaoglou et al. 



Wvithe & Loeb 



and sem i-analytic 



2003 ; Hatziminaoglou et al 
( see, e. g.. [cattaneo^^l| 
1200(1 ICavaliere & Vittori: 



2001 
2003 
1999; 



Enoki et al 



2000; 



2003 



Kauffmann & Haehnelt 

I Cattaned 1200 it ICav aliere & Vittorinil 120021 
Ivolonteri et all bo03t ISpringel et alJ l2005t ICattaneo et alJ 12005 ) 
hierarchical models have been proposed to understand the QSO 
phenomenon. Most of these models successfully reproduce the 
QSO luminosity function [LF] at high (i.e. z > 2) redshifts, when 
the QSO activity reaches its peak, and some of them also account 
for the existence of QSO at very high (i.e.. z ~ 6) redsh ift, hosting 
BHs as massive as 10 9 M S (Volon teri & Reesl 120051) . However, 
hierarchical models generally fail to match the QSO properties 
at lower redshifts since the halo merger rate, which is generally 
regarded as the mechanism that triggers the QSO activity, declines 
with time less rapidly than the observed QSO number density. Var- 
ious possible solutions have been proposed to solve this problem. 
For example it has been suggested that encounters within galaxy 
groups and clusters play a fundamental role in triggering the BH 
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accretion iCavaliere & Vittorinil boOO). that either the fraction of 
accreted ga s or the time-scale of acc retion depend on redshift 
lHaiman & MenouteOOollCattaneofcOOlh . or that the accretion rate 
is pro portional to the gas density of the host galaxy ICattaneo et alj 
2005). These solutions improve, with a different degree of success, 
the match to the observed LF at low redshift. The model proposed 
by Enoki et al. 2003 also allows to reproduce the QSO clustering, 
but is not guaranteed to match the correlation of Ferrarese (2002) 
between BH and halo mass. Another possibility, suggested by 
the outcome of hi gh resolution, hydrodynamical simulations 
I Hop kins etai]E o05). is that of assuming that for a large fraction 
of its lifetime the QSO is heavily obscured by surrounding gas and 
dust except for a short window time of ~ 10 7 yr during which it 
shines at its peak luminosity. 

In this work, we investigate to which extent 'standard' ana- 
lytic and semi-analytic hierarchical models, where the assembly of 
supermassive black holes is related to the merger history of DM 
haloes, can reproduce both the observed QSO LF and their cluster- 
ing at low redshifts. Fo r this purpose we imple ment the two orig- 
inal Mial^icjnodels ofJ^i&^^^e3 j2002) (WL02 hereafter) 
an d Wvit he & Loebl 12003^ (WL03 hereafter) and the semi-analytic 
one of Volonteri, Haardt & Madau (2003) (VHM hereafter), com- 
pare their predictions to the LF and clustering of the optical QSOs 
measured in various redshift catalogues, and discuss possible mod- 
ifications to the models that improve the match to the data. 

The outline of the paper is as follows. In Section 2 we present 
the observational datasets and the QSO properties relevant for our 
analysis. In Section 3 we present the WL02 and WL03 models, 
compare their predictions with observations and introduce some 
modification to the original models to better match the observed 
LF. A similar analysis is repeated in Section 4 for the semi-analytic 
model of VHM. Finally, in the last section we discuss our results 
and draw our main conclusions. 

Throughout this paper we assume a flat ACDM cosmological 
model with Hubble constant h = Hq/IQQ kms^'Mpc -1 =0.7 and a 
dominant contribution to the density parameter from the cosmolog- 
ical constant, Slj^ = 0.7. We adopt a CDM density power spectrum 
with primordial spectral index n = 1 and normalized by assuming 
o 8 =0.9. 



2 DATASETS 

The main dataset considered in this work is the QSO redshift cat- 
alogue by ICroom et alj 120041) [C04 hereafter] obtained by merg- 
ing the 2dF QSO Redshift Survey (2QZ), containing objects with 
an apparent bj magnitude 18.25 < bj < 20.85, with the 6dF QSO 
Redshift Survey (6QZ) of bright (16 < bj < 18.25) quasars. The 
full sample includes 23660 quasars in a wide redshift range (0.3 < 
z < 2.9) and spread over 721.6 deg 2 on the sky. The 2QZ/6QZ cat- 
alogue is affected by various types of incompleteness described in 
details by C04 that need to be accounted for in order to minimize 
systematic effects. For this purpose we have considered various 
subcatalogues, described below, characterized by a higher degree 
of completeness to allow a precise measurement of the QSO lumi- 
nosity function [LF] and clustering properties. 

A second dataset consisting of 5645 quasars with an apparent 
magnitude 18.0 < g < 21.85 extracte d from the recen t 2dF-SDSS 
LRG and QSO [2SLAQ] survey iRichards etail2 005. R05 here- 
after) has also been considered to provide us with a reliable esti- 
mate of the LF of fainter objects at z < 2. 1. 



2.1 The QSO Optical Luminosity Functions 

The 2QZ/6QZ catalogue has been used by C04 to compute the op- 
tical LF of a subsample of 15830 QSOs brighter than M b . = 22.5 
in the redshift range 0.4 < z < 2.1. The cut in absolute magni- 
tude guarantees a minimum spectroscopic sector completeness of 
at least 70 per cent, while redshift constraints ensure a photomet- 
ric completeness of 85 per cent. The LF has been evaluated into 
8Mhi = 0.5 bins in absolute magnitude using the 1 /V estimator of 
IPage & Carreral fcOOOl) into six equally spaced, independent red- 
shift bins. The filled dots in the six panels of Fig.Qshow the LF of 
C04 measured in the different redshift intervals indicated in each 
plot. 

In the same figure the thin solid line represents the best fit to 
the 2SLAQ LF determined by R05 in the same redshift intervals 
and bj magnitude bins. The original g magnitudes have been trans- 
formed into bj magnitudes by using the relation (g — bj) = —0.045 
of R05. Finally, bj magnitudes have been converted into B-band 
magni tudes by using the relation Mb = Mj,j + 0.07 o flPorciani et alj 
(2004) [PMN hereafter]. 



2.2 QSO vs. Dark Matter Clustering 

PMN have estimated the QSO two-point spatial correlation func- 
tion of ~ 14000 2QZ/6QZ quasars with redshift 0.8 < z < 2.1 in 
three different redshift intervals [0.8,1.3], [1.3,1.7] and [1.7,2.1]. 
The three subsamples with median redshifts z e ff = 1.06, 1.51, 1.89 
contain ~ 4300, ~ 4700 and ~ 4900 objects each. The more con- 
servative cut in redshift and the use of these three redshift intervals 
guarantee (i) a photometric completeness larger than 90 per cent, 
(ii) a similar number of quasars in each redshift bin and (iii) that 
each subsample covers a similar interval of cosmic time. 

PMN have measured the spatial two-point correlation func- 
tion, ^ obs (r ± ,7t), using all QSO pairs separated by 7t along the line 
of sight and by along the direction perpendicular to it and have 
computed the 'projected correlation function' S by integrating ^ obs 
along the n direction: 



S obs (rJ 



r± Jo 



(r±,n)dTZ . 



(1) 



Then, in order to quantify the QSO clustering with respect to the 
mass, they have evaluated the quasar-to-mass biasing function 

£. obs (r ^,Zmin < 2 < £max) 



(2) 



EDMf/X.Zeff) 

where z m ; n and z max represent the lower and upper limits of the red- 
shift interval and the projected mass correlation function, 3dm (n)> 
was estimated as in Peacock & Doddj ll99a) . PMN have found 
that the QSO biasing function in Eq.0 is almost independent of 
the projected separation r±_ and thus it is possible to characterize 
the QSO spatial correlation properties at Zeft using a single 'bias' 
parameter b(z e g). The results of the PMN analysis are presented 
in Fig. [2] where the filled dots show the value of i(z e ff) for the 
2dF/6dF QSOs at three different redshifts, together with their lo 
uncertainty. The QSO-to-mass bias parameter b(z e s) increases with 
redshift, in quantitative agreement with the results of a similar anal- 
ysis performed bv lCroom et al J 12005). 



3 ANALYTIC MODELS FOR THE QSO EVOLUTION 

In this section we review the two analytic models by WL02 and 
WL03, use them to predict the QSO LF and clustering at z < 2 and 
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Figure 1. The QSO luminosity function in B-band at six different redshifts: models vs. observations. The filled circles show the 2dF/6dF QSO luminosity 
function measured by C04 together with their la error bars. The thin solid line shows the best-fit to the 2SLAQ QSO luminosity function of R05. The short- 
dashed and the dotted lines show the WL02 model predictions obtained by setting (y = 4.71, E = 10~ 5 1 , td c ,o = 10 6 3 yr) (label: WL02A) and (y= 4.71, e = 
10~ 5 1 , tdcfl = 10 7 2 yr) (label: WL02B), respectively. The thick solid line shows the WL03 model predictions obtained for (y = 5.0, e = 10~ 5 ' 7 ), while the 
thin long-dashed line (label: WL03k) represents the same model with an exponential cut in the luminosity-mass relation as in eg. 1141 . with k = 0.77 and M^. ilo 
indicated in the plots. 



compare the results to the observations described in Sections 12.11 
and 12.21 Since we apply these models beyond the redshift range 
they were originally designed for, we discuss and introduce some 
modifications to improve the fit to the 2dF/6DF and 2SLAQ LFs 
and clustering. 

3.1 The WL02 model 

The WL02 model describes the QSO evolution within the theoret- 
ical framework of hierarchical structure formation. However, in- 



stead of associating the QSO activity to the formation of DM viri- 
alized haloes, as originally proposed by Haiman & Loeb 1 1998), it 
assumes that the QSO phenomenon is triggered by halo-halo merg- 
ers. The model also assumes that the mass of the BH powering the 
QSO, Mbh, is a fraction, £, of the host halo mass Mj, a i and that, 
after a merging event, the QSO shines at the Eddington luminosity, 
^Edd> w ith an universal light curve, /(?). The S-band QSO lumi- 
nosity can be related to My, and Mh a i through /(f): 

MO = M bh/(0 = eMhalo/O) for M halo > M min , (3) 
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Figure 2. The mean QSO-to-mass biasing parameter, b(Zetf) estimated at 
three effective redshifts z = 1.06, z = 1.51 and z = 1.89: models vs. ob- 
servations. The filled circles show the mean biasing of 2dF/6dF quasars 
measured by PMN with the associated 10 uncertainties. The short-dashed 
line shows the theoretical predictions of the WL02 model with (y = 
4.71, e = 10 _5 - 1 , fjc.o = 10 6 ' 3 yr) (label: WL02A) and with (y=4.71, E = 
10~ 51 , t dc , = 10 7 2 yr) (label: WL02B). The amount and evolution of QSO 
biasing of these two models are exactly the same. The thick solid line shows 
the WL03 model predictions obtained for (y = 5.0, £o = 10 -5 ' 7 ) and the 
thin long-dashed line (label: WL03k) represents the same model with an 
exponential cut in the luminosity-mass relation as in eq. fl4l. with k = 0.77 
and M|* ato indicated in the plot. 



where M m ; n ~ 10 8 [( 1 + z)/10]~ 3 / 2 M© is the minimum halo mass 
inside which a BH can form. The number density of active QSOs 
then can be obtained by multiplying the number of haloes with 
mass between AMh a i and AMh a i + dAM b . d \ that: accrete onto a 
halo of mass Mh a i — AMj, a i per unit time by the number density 
of haloes in the same mass range: 



/(M halo ,AM hal0 ) = 



dN(M, 
dM~ 



IM=M M 
d 2 N a _ 



-AM hi 



dAM hdlo dt lM=M hal0 -AM hal , 



(4) 



the halo mass function dN j^' z 



The quantity J^"'" 8 ^ represents the me rging rate predi cted by the 
Extended Press-Schechter theory [EPS] Lacev & Cole 1 1993) and 

J 19991) . 



is that of Sheth & Torme 
A further relation between My, and Mh a i can be imposed by 
assuming that the scaling relation between the BH mass and the 
circula r velocity of the host DM halo, v c , observed by Ferrarese 
(2002) in the local universe, holds at all redshifts: 



M bh « v l = (159.4)1 



M ha lo 



W u h- l M Q 

Qm(0) 



7/3 

A c \r/6 

18712/ 



(l+z) T/2 , (5) 



where the second equality follows from the relation between v c 
and M halo {Bark ana & Loebl200ll) in which A c (z) = 187t 2 + %ld - 
39d 2 , d = £2 m (z) — 1 and £2 m (z) represents the matter density pa- 



rameter at a given redshift z. From eqs.J3J and J5j: 

°Uo 12 M ; \£i m (z) 18**/ { ' ' W 
Finally, we assume that the QSO luminosity curve is given by 
a simple step function 

iEdd,Sg/'AM halo 



-tdcfl - 



(7) 



where f^ c o <C// 1 (z) is the time of QSO duty cycle at z = 0. 
Eqs.{5J, {5J and Q allow us to compute the QSO LF: 



*(L B ,z) = f" dM hh [ 



0.5eM Wo 



dz 



EM m 
,dN hh \ 



dAM bb 

d ./V mersn 



dM I M=Afbh — AMbh dAMfrfrdt lM=Mbh— AM b h 

dt' 

x—S[L B -M hh f(t z -t')] 
dz 



(8) 



that, once integrated over M D h, has an analytic expression that de- 
pends on the three free parameters /j c q, £q and y: 



*(**.z) 



./o 



0.5M halo 3 
rfAM halo - 



f rfc n AM ha 



halo 



'y£5.7x 10 3 M hal0 

x / (M ha io, AM halo ) , (9) 

where M ha , = L Edd , B /(5.7 x lO 3 eL 0iB )M e . 

The connection between QSO luminosity and halo mass 
in our model, Ls(M b . d i ), and the existence of analytic models 
that describe the spatial and clustering of DM haloes (see, e.g., 
iMo & Whitdll996l:ISheth & Tormenlll999h allow us to investigate 
the spatial correlation properties for QSOs. In particular, we can 
compute the QSO-to-mass biasing parameter, b(z), averaged over 
all QSO luminosities: 



b(z) 



b(L B (M ha]o ) , z)*(L fl) z)dL B 



<S>{L B ,z)dL B 



(10) 



where S>(L#,z) is the QSO LF in Eq.{5), L B m [ n is the luminosity 
of the faintest object in the sample. The quantity b(L B (M bli i ),z) 
represents the biasing parameter of a halo of mass Mh a i hosting a 
QSO of lu minosity L B at the redshift z that has been computed by 
ISheth&Tormenljl999t) : 



&(MM halo ), Z ) = l + -i-[^-l| 



+ 



(11) 



5 C (0) Vl + [^8 e (z)/o M ] 2 ^ 
where a = 0.707, p = 0.3, 8 c (z) is the critical threshold on the lin- 
ear overdensity for spherical collapse at redshift z and aj^ is the rms 
linear density variance smoothed with a 'top-hat' filter correspond- 
ing to the mass M. Eq. jlOt provides us with an analytic expression 
for b(z). It assumes a univocal relation between quasar luminosity 
and halo mass, or, in other words, that the probability of finding a 
quasar of a given luminosity L B in a halo of mass Mh a i depends on 
Mhalo only- 

3.2 WL02 model vs. observations 

Let us now compare both the QSO LF and the biasing parame- 
ter predicted by the WL02 model to the various datasets avail- 
able. Given the background cosmology, the model predictions are 
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fully specified by a set of three parameters: (y, £o, tdc.o)- Here 
we explore two separate cases. The first one, which is slightly 
different from the one considered in the original WL02 paper 
(and labeled WL02A in Figs. □ and |2}, uses the parameters (y = 
4.71, £ = 10" 51 , t dcfi = 10 63 yr). The firs t two chosen val - 
ues represent the best fit to the observations of iFerraresd 120021) . 
while the latter, which corresponds to a value of f^ c (z = 3) = 10 6 ' 9 
yr at z ~ 3, is fully consistent with the lifetime of bright QSOs, 
t c / c (z — 3) = 10 7 yr, inf e rred f rom the sample of Lyman-break 
galaxies of ISteidel etall <2002h . The predicted QSO LFs, plot- 
ted as short-dashed line in Fig.Q fail to match the observed LF 
both in the bright and the faint ends. This result is similar to 
the one originally obtained by WL02 using the set of parameters 
(y= 5.0, e = 1(T 54 , tacfi = 10 6 ' 3 y). that constitutes their best 
fit to the data. Since the overall amplitude of the model LF (Eq. 
|9j linearly depends on y _1 , Eg 1 and tj c Q, it is quite straightfor- 
ward to boost up the model LF to match the number density of 
the observed QSOs. For example, fixing the values of y = 4.71, 
£o = 10 , and leaving fycO as a free parameter, we find a best fit 
for td c o = 10 7,2 yr. The resulting model, labeled WL02B, is shown 
in Fig.Qas a dotted line. This duty-cycle is very large and some ad 
hoc modifications to the WL02 model would be req uired to satisfy 
the high-redshift constraints of Steid el et alj (2002). We will dis- 
cuss physically motivated modifications to the WL02 model later 
on in the framework of the WL03 model. 

The biasing functions predicted by the WL02A and WL02B 
models are shown in Fig. |2|at three different redshifts and com- 
pared to the PMN data. The line-styles are the same as in Fig.Q 
The two models predict the same amount and evolution of QSO 
biasing. In both cases at z = 1.06, QSOs are mildly biased with 
respect to the underlying mass, marginally consistent with obser- 
vations, while their clustering at z = 1.89 is significantly less than 
observed. 



3.3 The WL03 model 

WL03 modified the original WL02 model by prescribing a 
self-regulated accretion mechanism of SBHs following major 
mergers betwee n haloe s of different masses, as proposed by 
Kauff mann & Haehneltl fcOOOl) . A self-regulated accretion mech- 
anism is the natural outcome of the production of powerful gas 
winds that interrupt the infall of gas on the BH after halo mergers. 
Self-regulation takes place when the energy in the outflow equals 
the gr avitational binding energy in a dynamical time ISilk & Reesl 
1998). Assuming that the gas is located in a disk with characteristic 
radius ~ 0.035r v ; r , the dynamical time /dyn is given by 



dyn 



= 0.035- 



3.64 x 10 7 ( a m (Q) A c 



-1/2 



Q m (z) 18ti 2 

x(l+zr 3/ V, (12) 

and represents the quasar duty cycle: t dc = f dyn . The major merger 
condition has been introd uced to gu arantee that the dynamical 
friction time-scale iBinnev & TremaineMl987l) for the satellite is 
shorter than a Hubble time. For this reason WL03 only consid- 
ered mergers between haloes with a mass ratio P = AMh a i /Mh a i > 
0.25. As a consequence, the model LF can be expressed as 



^dyn 



fO.SMralo 3 

<S>(L B ,z) = dAM h . dio — 

J0.25Mmo ye 5.7 x 10- 3 



x/(M halo ,AM halo ) . (13) 



It only depends on two free parameters: £o and y. The model biasing 
function has been computed using Eq. <10> . 



3.4 WL03 model vs. observations 

The first model we have explored is the one by WL03 correspond- 
ing to the choice of parameters (y = 5.0,£n = 10~ 5 7 ), still con- 
sistent with the observational data of lFerraresel 120021) . The model 
LF (labeled WL03 and plotted with a thick solid line in Fig.0 is 
very similar to that of WL02B, but has the advantage of having a 
physically, well motivated QSO duty-cycle. The WL03 LF matches 
the observed one at low luminosities, but it overpredicts the num- 
ber density of bright QSOs. This discrepancy is more evident at 
low redshifts and can be accounted for by modifying the WL03 
model. One possibility is to allow for a major merger threshold P 
that depends on Mj, a i . As we have checked, it is indeed possible 
to find some suitable function P{M^. i \ Q ) monotonically increasing 
with M}, a i that allows to match both the faint and bright ends of 
the observed QSO LF. While this is somewhat an ad hoc solu- 
tion, a more physically plausible modification has been proposed 
by WL03 and consists of assuming that accretion onto BH is ham- 
pered by the high temperature of the gas within group/cluster-size 
haloes. WL03 proposed that accretion onto the central BH should 
be prevented within haloes of masses larger than 1O 13,5 M0, corre- 
sponding to a L ~ 2 x IO^Lq b at z — 1. We have made a simi- 
lar assumption and proposed that the accretion efficiency decreases 
above a given critical halo mass resulting in a modified relation 
between Lqso and Mj, a i : 

L QS o=L(\-exp(-(L*/L) k )) . (14) 

In the previous equation L ===== 5.7 x 10 3 £(Mh a i o /M Q ) is the B-band 
Eddington luminosity of the original WL03 model and L* is the 
Eddington luminosity of a halo with critical mass M^ alo . Both k 
and Mj^ alo can be treated as free parameters. The thin dashed line in 
Fig. Q shows the effect of a keeping k = 0.77 while leaving M^ fl/o 
as a free parameter, whose value is indicated in the plots. The LF 
predicted by the model (labeled WL03k) provides a good match 
to the observed one at all redshifts, except perhaps for very bright 
objects. The resulting values of M^ alo range between 10 129 and 
10 13 ' 4 M a which constitute plausible values, close to that proposed 
by WL03. Leaving also k as a free parameter does not improve 
significantly the agreement with the observational data and a best 
fit value close to k = 0.77 is found at all redshifts. We have also tried 
to use the model WL03 without the restriction to major mergers, but 
it has resulted inconsistent with the data by about several orders of 
magnitude and we have decided not to show it. 



4 SEMI-ANALYTIC MODELS FOR THE 

HIERARCHICAL EVOLUTION OF QUASARS 

Semi-analytic approaches to the evolution of QSOs in a hierarchical 
scenario also assume some relation between QSOs and DM haloes 
properties. However, they differ from analytic modeling since the 
evolution of DM haloes and of the QSOs within them are treated 
separately. The merging history of DM haloes is described by the 
EPS formalism. Phenomenological relations are used to model the 
physical processes leading to the quasar evolution. It is therefore 
possible to adopt a more detailed description of the physics involv- 
ing the baryonic component of cosmic structures, including BHs. 
A second advantage of the semi-analytic approach is the flexibil- 
ity of the scheme, so that different astrophysical prescriptions can 
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be tested within the same framework of cosmic evolution of DM 
haloes. 

4.1 The VHM model 

In this work we focus on the semi-analytic model developed by 
VHM that describes the hierarchical assembly, evoluti on and 
dynam i cs of the BHs powering QSO s I Volonteri, Haardt & Madau 
| 2003l IVolonteri Madau & Haardtl 120031 iMadau et alj 120041 
Ivolonteri & Reed I2005I) . Like WL02 and WL03, VHM also 
assume that: (i) the observed correlation between BH masses and 
circular velocity (Ferrarese 2002) justifies the assumption of a 
link between QSO activity and haloes' properties and constitutes 
a constraint to the semi-analytic model at z = 0, and that (ii) the 
cosmological evolution of DM haloes is well described by the EPS 
theory. Moreover, like in the WL03 case, they further assume that 
(Hi) the QSO activity is triggered by major mergers. 

The numerical implementation of the semi-analytic model 
consists of a two-step procedure. The first step is aimed at con- 
structing a set of halo merging histories using the EPS theory. In 
the EPS formalism, when one takes a small step 8z back in time, 
the number of progenitors a parent halo of mass Mq at z = Zo frag- 
ments into is lLacev & Colell993l) : 

,2 



dN 1 Mq 1 d8 c da 



M 



Sz, 



(15) 



where S = <3 2 M (z) - a^ (zq), G 2 m {z) and a 2 Mg (zq) are the linear the- 
ory rms density fluctuations smoothed with a 'top-hat' filter of mass 
M and Mq at redshifts z and zo, respectively, and 8 c (z) is the crit- 
ical thresholds on the linear overdensity for spherical collapse at 
redshift z. Integrating this function over the range < M < Mq 
gives unity: all the mass of Mq was in smaller subclumps at an ear- 
lier epoch z > zo- From Eq, J15> we can compute a fragmentation 
probability that, via rejection methods, can be used to construct a 
binary merger tree. Implementing a successful Monte Carlo proce- 
dure, however, requir es the use of two different numerical approx- 
imations ISomerville & Kolatdf l 999 ) . First of all, since in a CDM 
cosmology the number of haloes diverges as the mass goes to zero, 
it is necessary to introduce a cut-off mass, M res , that marks the tran- 
sition from resolved progenitors (having M > M res ) to the accreted 
mass that accounts for the cumulative contribution of all mass ac- 
creted from unresolved haloes. Secondly, the time-step Sz has to be 
small enough to guarantee a small mean number of fragments N„ 
in the range M res < M < Mq /2, to avoid multiple fragmentation. 

Once the appropriate choices for 8z and M res are made, the 
binary tree i s constructed by u sing a Monte Carlo procedure simi- 
lar to that of ICole et all fcOOOl) . described in VHM. We have taken 
M re s = 10~ 3 M ( ) at z = decreasing with redshift as (1 + z) 3,5 . Fi- 
nally, for each Monte Carlo realization we have used 820 time- 
steps logarithmically spaced between z = and z = 20. As shown 
by VHM, with this parameter choice our merger tree algorithm not 
only conserves the mass, but also reproduces the EPS conditional 
mass function at all redshifts. 

In the second step of the procedure, we implement a set of 
analytic prescriptions and follow the accretion history of massive 
BHs within their host haloes to model the QSO activity. The VHM 
model assumes that the seed BHs formed with masses of 150M : ;, 
(note that, as shown by VHM, the final results are not very sensitive 
to this choice) following the collapse of the very rare Pop III stars, 
in minihaloes forming at z = 20 from the density peaks above a 
3.5o" threshold. In the assumed ACDM cosmology this corresponds 
to minihaloes with mass ~ 1.6 x lO 7 M .Then we assume that the 
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Figure 3. Fraction of haloes hosting at least a nuclear massive BH vs. red- 
shift. The long-dashed curve shows the occupation fraction computed by 
weighting over all branches of the merger trees. The occupation fraction 
increases with increasing halo mass: A/h a i > 1O IO M0 (dot-dashed curve), 
^halo > 1O U M (short-dashed curve) andM|, a i > 10 12 Mr., (solid curve). 



quasar activity is triggered only by major mergers with P > 0.1, a 
threshold lower than P > 0.25 adopted by WL03. 

Two main features differentiate the VHM model from the 
WL03 one. First, the VHM model is naturally biased, as the BH 
seeds are associated with high-density peaks in the fluctuations 
field. Second, VHM take into account the dynamical evolution of 
BHs, including strong gravitational interactions such as the gravi- 
tational rocket effect {Merri tt et al 12004b IVolonteri & Pernal 2005). 
Such dynamical interactions can possibly eject BHs at high veloc- 
ities from the centre of haloes. The net effect is to contribute to 
selecting massive haloes (i.e. those with a large escape velocity) 
as BH hosts. In all three scenarios we consider (see below), we 
have included a treat ment of t he 'gravitational ro cket' effect fol- 
lowine lFavata et alj 120041) and lMerritt et alj <2004 (upper limit to 
the rec oil velocity). More details can be found in Volon teri & Pernal 
(2005). Fig. [3] shows the occupation fraction, i. e. the fraction of 
haloes hosting nuclear BHs for haloes of different masses. The oc- 
cupation fraction of large haloes (Mj, a i > 1O 12 M0) is of order unity 
at all times, while smaller haloes have a large probability of being 
deprived of their central BH. 

Following a major merger, the BH at the centre of the mas- 
sive progenitor can grow in mass in two ways: (i) after a dynam- 
ical friction time when a bound BH binary system form s at the 
centre of the halo, hardens via three body interactions lOuinlanl 
1996; Milosavlievic & Merritt 2001) and then rapidly coalesces 
through the emission of gravitational waves (Peters 1964); (ii) af- 
ter a dynamical free-fall time when a signi ficant fraction of the 
gas falls to the centre of the merged system I Spring e! et aljEo05l 
iDi Matteo et all2 005) and is accreted on the BH at an appropriate 
rate. Yet, as shown by VHM, the first mechanism contributes little 
to the BH mass accretion and will be neglected in this work. To im- 
plement the second mechanism we need to specify the prescription 
for the mass accretion and its rate. We have explored three different 
scenarios. The first two assume that BHs start accreting mass at the 
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Figure 4. The QSO luminosity function in B-band at six different redshifts: models vs. observations. The filled dots show the 2dF/6dF QSO luminosity 
function measured by C04 together with their la error bars. The thin solid line shows the best fit to the 2SLAQ QSO luminosity function of R05. The open 
squares refer to the VHM model predictions when a fraction a = 7 x 10~ 6 of accreted mass of the merged system total mass is assumed. The open triangles 
show the VHM model predictions when a relation between the accreted mass and the circular velocity of the host halo is assumed. 



Eddington rate after about one dynamical free-fall timescale from 
the merger. Accretion lasts until a mass AM accr has been added to 
the BH, but, as in VHM, we have inhibited gas accretion in all 
haloes with v c > 600 kms -1 . 

In the first of the two scenarios, labeled El, the accreted 
mass is proportional to the mass of the available gas and hence to 
the total mass of the massive progenitor: AM accl - = aM^. d \ . Here, 
oc = 7 x 10~ 6 guarantees the normalization of the — a g re- 
lation at z = 0, where g g is the velocity dispersion of the host 



galaxy iTremaine et alj|2002|). scaling with the circular velocity 

No feedback is ex- 



galaxy I Iremaine et al. 2002), scaling wttl 
of the halo as suggested bv iFerraresel 1200! 



plicitly included in the semi-analytic scheme. As Mj, a i °= vj?, the 
slo pe of the Mm, — <3g rela tion is flatter than the observed one (but 
see Wvithe & Loeb 2005a). This scenario is similar in spirit to the 
WL03 model, and is meant to compare the clustering properties of 
quasa rs at low redshift to th at of their higher redshift counterparts. 
Adelberger & Steidel 1 2005) indeed find that the clustering of ac- 
tive nuclei at 2 < z < 3 points to a — M^ a \ relation which is 
independent of redshift. 

The second scheme [E2] assumes a scaling relation between 
the accreted mass and the circular velocity of the host halo, 
AM accl - = kxv^, which is normalized a-posteriori to reproduce the 
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Figure 5. The QSO luminosity function in B-band at six different redshifts: models vs. observations. The filled dots show the 2dF/6dF QSO luminosity 
function measured by C04 together with their la error bars. The thin solid line shows the best fit to the 2SLAQ QSO luminosity function of R05. The open 
pentagons refer to the VHM model predictions including super-critical accretion in high-z haloes. The asterisks show a model in which the accreted mass 
scales with the circular velocity of the host halo, and accretion is suppressed in haloes with Mi, a i > ID 1 x < \t . . 



observed relation between M\,^ and v c at z = jFerraresefe ot^l. 
VHM showed that this scenario overestimates the optical LF at 
z < 1. Here, we assume a linear dependence of k on redshift, as 
k = k(z) — 0.15(1 + z) +0.05, to account for the decrease of the 
gas available to fuel BHs. The above relation is totally empirical 
and was determined iteratively in order to get a better fit to the LF, 
without underestimating the local SBH mass density. 

The last prescription for mass accretion, labeled B, assumes an 
early stage of super-critical accretion during which the central BH 
accretes mass at a rate that can be estimated by the Bondi-Hoyle 
formula !Bond i&HovIelll944l) . This model applies to metal-free 



haloes, therefore we assume that by z = 12 the interstellar medium 
has been enriched, and we inhibit super-critical accretion rates. 
When the super-critical phase ends, accretion proceeds in subse- 
quent episod es as in model E2. This p ossibility has been recently 
advocated by I Volonteri & Re"esl (2005) to reconcile a hierarchical 
evolution with the existence of QSO at z ~ 6, hosting SBHs with 
masses ~ 10 9 Mq. 

The end product of our semi-analytic models is a set of merg- 
ing and accretion histories for 220 parent haloes with masses in the 
range (1.43 x 1O H M , 1O 15 M ). When active, i.e. during the pe- 
riod of mass accretion, the QSO shines with a B-band luminosity of 
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(Lb/Lq) = Mbh x 1O 3 46 M , obtained under the assumptions that 
the rest mass is converted to radiation with a 10 per cent efficiency 
and that only a fraction /b = 0.08 of the bolometric power is ra- 
diated in the blue band. Finally quasars in the model outputs are 
selected according to the 2dF/6dF criteria. 

The model QSO LF at different redshifts has been computed 
by evaluating the number density of active QSOs in each luminos- 
ity bin in redshift intervals centred on an effective redshift z en -. In 
practice we have counted the number of active QSOs in the redshift 
and luminosity bins in all merger trees, each of them weighted by 
the number densi ty of their parent haloes at z = [evaluated using 
the S heth & Tormenl ^999) formula]; the result has been normal- 
ized using the number of time-steps in the redshift interval and the 
number of merger trees considered. Associate uncertainties have 
been computed by assuming Poisson statistics. 

The biasing function b(z) at the three redshifts considered by 
PMN has been estimated by using the following equation 



Jo 



/ «P(A4alo(£fl > L - t n 
JO 



, (16) 



,z)dM h 



where *P(M,z) is the mass function of the haloes hosting 
QSOs with luminosities larger than the selection thresholds of 
PMN (L B , min /L ,B = {1.5 x 10 n ,3.9 x 10 n ,6.6 x 10 11 }), and 
b(Mb : , ]n,z) is the bias function of haloes computed following 
ISheth&Torrner]<1999h . 

Model predictions also allow a straightforward evaluation of 
the mean halo occupation number, i.e. the average number of ac- 
tive QSOs hosted in haloes with mass between Mj, a i and Mj, a i + 
rfMhai : 



NQSo{Mhaio) = 



dN/dMi 



(17) 



halo 



where ^(M^aio) is the mass function of haloes hosting ac- 
tive quasars in the h alo merger trees and dN/dM\ rd \ is the 
ISheth&Torrner]< 19991) halo mass function. Model uncertainties 
have been evaluated from Poisson errors associated to ^(Mhaio). 



4.1.1 VHM model vs. observations 

In Figs. [4] and [5] we compare our model LF with the observations. 
The 6dF/2dF and the 2SLAQ QSO LFs are plotted using the same 
symbols as in Fig.Q The LFs predicted by our semi-analytic mod- 
els are plotted with the three different symbols indicated in the 
plot. At low redshifts all models reproduce the faint end of the LF 
fairly well. However, they systematically overpredict the number 
of bright QSOs: this indicates that having inhibited gas accretion 
in haloes with v c > 600 kms -1 has little impact on our results. 
Indeed, this circular velocity is significantly larger than that as- 
sociated to the mass threshold A/^ alo ~ 10 13,5 Mq adopted in the 
WL03K model. Imposing a similar mass threshold in the VHM 
model (model E2 cut ), in place of the v c > 600 kms -1 cut-off, we 
find a much better agreement between the model predictions and 
the observed LF at all z < 1.5. Finally, we note that model El fails 
to produce QSOs brighter than lO 12 /^^. 

As shown in Fig. [6] all semi-analytic models explored match 
the observed biasing function out to z ~ 2. We notice that the three 
semi-analytic models ex plored predict a very sim ilar biasing func- 
tion. As pointed out by IWvithe & Loef] 1200 5a). the evolution of 
clustering is slightly faster when the BH mass scales with the halo 
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Figure 7. The mean halo occupation number of QSO at three effective 
redshifts z = 1.06, z = 1.51 and z = 1.89 (from top to bottom). Histograms 
represent model predictions, while the dashed line shows the mean halo 
occupation number proposed by PMN and consistent with 2dF/6dF QSO 
data. Left panels: model El. Central panels: model E2. Right panels: model 
B. 



mass (model El) rather than the circular velocity (model E2), al- 
though the difference is of little significance given the large scatter 
in the model predictions. 

The reason of the difference can be understood by looking at 
the histograms plotted in Fig.0representing the mean halo occu- 
pation number, A'Qso(^halo)' °f tne various models. At z = 1.89 
model El exhibits a steeper dependence on halo mass compared 
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Figure 8. Accreted mass as a function of redshift for different halo masses. 
From top to bottom, curves are for 2.5a, 2a and 1.5a peak haloes. Dashed 
curves: model El. Solid curves: model E2. The horizontal line shows the 
BH mass corresponding to LB, m i n at z = 1.89 of PMN, assuming Eddington 
accretion rate. 



to the other cases meaning that QSOs are preferentially found in 
massive haloes with a high degree of spatial clustering. This re- 
sult derives from the fact that the accretion scheme of model El is 
more efficient at very high redshift (z > 6) than that of model E2, 
due to the different scaling of AM accl - with redshift. Fig. [^exempli- 
fies this effect. On one hand, AM accr is a steeper function of redshift 
in model E2, implying that massive BHs accrete more mass in ev- 
ery accretion episode, thus leading to a longer duty-cycle, and a 
larger occupation number, in general. On the other hand, at redshift 
Z = 1.89, only quasars above L B m [ n are selected. In model E2 the 
2<T-peak haloes contain BHs massive enough to be above threshold, 
while for model El a slightly more massive halo is needed. Con- 
sequently this enhances the bias. At lower redshifts this effect be- 
comes prog ressively less importan t . This result interestingly agrees 
with that of Adelbereer & Steidel 1 2005): at z > 2 SMBH masses 
correlate with the halo mass instead of with velocity dispersion, or 
circular velocity. We can speculate that there might be a transition 
in the interaction between SMBHs and their hosts, which switches 
on at z — 2: at higher redshifts the SMBH mass scales with the halo 
mass, at lower redshifts with its velocity. As shown in Fig.[S| the ac- 
creted mass in a given episode is larger in model El than in model 
E2 for haloes representing density peaks below 2a at z > 3. Model 
El, therefore, implies an earlier growth for SMBHs. 

On the contrary, model B, in which massive BHs accrete mass 
with very high efficiency at high redshifts preferentially populate 
smaller haloes. Massive BHs in more massive haloes have already 
grown to masses close to the — a g threshold. At z = 1 .06 all 
models, especially El, predicts a number of QSOs in haloes which 
is systematically larger than the one inferred from PMN, especially 
in high-mass haloes. This discrepancy, which is marginally signifi- 
cant, considering the errors estimated by PMN, reflects the fact that 
semi-analytic models overestimate the optical LF of bright quasars 
at low redshifts, as shown in the top panels of Figs.|4]and[5] On the 
other hand, at high redshifts all models predict a halo occupation 



number that is systematically smaller than the PMN one, which 
again reflects the fact that our model LFs slightly underestimate 
the observed one at high redshifts (Figs. [4] and [5] bottom panels). 
This effect is particularly evident for model El, which, as we have 
noticed before, predicts no QSOs with L > 10 12 Lq b at z > 1.25. 



5 DISCUSSION AND CONCLUSIONS 

In this work, we have carried out for the first time a systematic 
study aimed at assessing how well, within the framework of hier- 
archical build up of supermassive black holes and massive halos 
hosts, analytic and semi-analytic techniques perform in reproduc- 
ing the QSO luminosity function and spatial distribution at very 
low redsh ifts, bound to reproduce the relation between My, and 
Mj, a io l Ferrarese 2002). We have found that only minor, physically- 
motivated modifications to standard semi-analytic techniques are 
required to match observations at redshifts as small as 0.5. More 
profound modifications seem to be required for a successful mod- 
eling of the very local QSO population. 

Previous studies, like those performed by WL02, WL03 and 
VHM, have shown that analytic and semi-analytic models well re- 
produce the observed QSO LF at high redshifts. At low redshifts, 
however, these methods systematically overpredict the number den- 
sity of bright objects; a mismatch that becomes increasingly large 
when decreasing the redshift. To improve the fit to the data we have 
proposed two simple modifications to the original WL02, WL03 
and VHM models. The first possibility is to assume that the ma- 
jor merger threshold required to trigger the QSO activity depends 
on the mass of the hosting halo. This reduces the number density 
of bright QSOs that can only be activated by rare merging events 
between large haloes of similar masses. A second possible mod- 
ification is physically motivated by the fact that gas accretion is 
expected to be inefficient within large haloes due to the high tem- 
perature of the baryons: this can inhibit accretion in haloes larger 
than ~ 10 13,5 Mq, hence reducing the abundance of bright objects. 
Both prescriptions significantly improve the fit to the bright end 
of the QSO LF, especially when implemented within the analytic 
framework. Yet, significant discrepancies still remain, especially at 
very low redshifts, that could be possibly eliminated by including 
two more factors that are missing from the models. 

First, when comparing our model predictions to the observed 
6-band QSO LF, we have implicitly ignored the presence of a sub- 
stantial population of (optically) obscured, luminous AGN at low 
redshifts, who se existence, instead, is sugg ested by Chandra re- 
sults (see, e.g., iBarser et alJ200lURosati et all2 002 ) . In our model- 
ing we have not included any correction for type-II qu asars, whic h 
appear to compris e between 30 per cent lLa F ranca et alj|2005|). 
and ~ 80 per cent (B rown et all2005llFranceschini et all20 05) of 
the quasar population at 0.5 < z < 2. However, our model suc- 
cessfully reproduces the QSO LF in the hard X-ray band. Fig. [5] 
shows a comparison between the model predictions and the [2-10 
keV] LF obtained bv lUeda et alj 1200 3l) . The conversion from the 
hard-X band luminosity into a blue-band estim ate ha s been done 
by assuming the scaling proposed by iMarconi et alj l2004) . We 
also compared the integrated mass density as a function of redshift 
to the redshift dependent integral of luminosity density of quasars 
iMarconi et al]|2004 . including the X-ray LF lUeda et alJl2003l) 
and a correction of a facto r 2 t o account for missing Compton- 
thick quasars I Br own et all2 005). We found a very good agreement 
over the whole redshift range probed by the LFs (0 < z < 3). We 
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Figure 9. Comparison of the model predictions to the observed optical QSO LF (filled circles: C04; thin line: R05) and the hard X-ray LF <Uedaetall2 ()()3). 
converted into B-band (see text for more details). Symbols are as in Figs.l4landl5l 



therefore ascribe the discrepancy between model predictions and 
observations at optical wavelengths mainly to selection effects. 



Second. lMerloni et al] (2003) and lMerlonil 12004) have shown 
that low-redshift AGN are probably accreting inefficiently, i.e. both 
at an accretion rate much smaller than the Eddington rate and with 
a low radiative efficiency. Evidence for a BH powering mechanism 
less efficient at low redshift is also provided by the fact that local 
bright QSOs seem to be hosted in early-type galaxies that show no 
sign of recent merging events like disturbed morphology or recent 
star formation episodes (s ee, e.g.. lGrogin et all2005tlGrazian et all 
12004 iDunlop et alj f2003 , and references therein). These consid- 
erations suggest that a successful model for describing the evo- 
lution of QSO luminosity should include both a prescription for 



a frequency-dependent galaxy obscuration and a more sophisti- 
cated mechanism for the QSO activity in which the BH accretion 
rate and the QSO duty-cycle might depend on halo masses and 
merger parameters, as suggested by recent numerical experiments 
iDi Matteo et all2005t) . 

Both analytic and semi-analytic models predict a moderate de- 
gree of QSO clustering at low redshift, consistent with the obser- 
vations. At higher redshifts (z ~ 2) the QSO biasing predicted by 
the analytic models appears to be significantly smaller than that of 
2QZ/6QZ quasars. The only way for reproducing the observed de- 
gree of clustering is to increase th e normalization constant e in the 
Mbh ~ v e relation (Wvithe & Loeb 2005a), which, however, would 
overpredict the QSO number density in the local universe. 
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On the contrary, the redshift evolution predicted by semi- 
analytic schemes is significantly faster and matches observations 
out to z ~ 2. The high degree of clustering predicted by the 
semi-analytic VHM model does not derive from having placed 
the first seed BHs in correspondence of high-o overdensity peaks 
of the mass density field. The adopted threshold, 3.5o at z = 
20, in fact corresponds to having at least one seed massive BH 
in haloes with M^\ ~ 10 M® at z = 0, which host massive 
BHs with masses well below that sampled by the optical LF of 
quasars . Moreover, dynamical effects such as the gravitational 
rocket IVolonteri & Pern Jl2005h can possibly eject BHs and thus 
lower the occupation fraction only in haloes with Mh a i < 1O 12 M 
that host BHs too faint to be included in the range probed by the op- 
tical LF. Placing seed BHs in correspondence of even higher peaks 
would certainly increase the biasing of the QSOs without modi- 
fying their luminosity function at z > 0.5, provided that the ma- 
jor merger threshold is changed accordingly (VHM). In this case, 
however, it would be difficult to explain the presence of SBHs in 
galaxies like the Milky Way or smaller, and, in general, AGN har- 
boured in dwarf galaxies < Barth et al. 2005), which anyway are not 
sampled by the quasar LF at z > 0.4. The large values of bias in the 
semi-analytic models have a different explanation: it derives from 
the lack of a deterministic relation between the DM halo masses 
and the QSO luminosities at a given time. Indeed, a finite time is 
required to accrete a mass AM accr to the central BH. During the 
accretion phase, the BH is smaller than predicted by the simple 
scaling relations with halo masses (i.e. — or M^h — -Mhalo)- 
This means that, on average, the hosting halo of a quasar of a given 
luminosity is larger in the semi-analytic scheme than in the ana- 
lytic models, the masses being the same only at the very end of the 
accretion episode. Consequently, the bias is enhanced in the semi- 
analytic model even if the LF looks similar. 

In our analysis we have shown that simple models in which 
QSO activity is triggered by halo mergers within the framework 
of hierarchical build-up of cosmic structures can quantitatively de- 
scribe the observed evolution of the quasar number counts and lu- 
minosity at all but very low redshifts, provided that some mech- 
anisms are advocated to inhibit accretion within massive haloes 
hosting bright quasars. Semi-analytic schemes, that naturally ac- 
count for the finite accretion time, are also capable of reproducing 
the observed QSO clustering at z < 2. These results suggest that 
the QSO evolution can probably be explained within the hierar- 
chical scenario of structure formation in which the QSO activity 
mainly depends on the masses of their host haloes. More realistic 
models, however, should also account for the various halo proper- 
ties including those of the baryons and for possible environmental 
effects. One possibility, which is most easily implemented within a 
semi-analytic scheme, is to include absorption effects and modify 
the accretion scheme accordingly using the outcomes of numerical 
hydrodynamical experiments. With this respect, we plan to modify 
the simple accretion schem e used in t he VHM model following the 
recent results o f lHopkinsetai]<2005h who found that active QSOs 
are heavily obscured for most of their time but during their peak 
of activity when the surrounding dust is blown away and the QSOs 
shine prominently for a short period (~ 10 7 yr). 

It is worth stressing that the hierarchical models for QSO evo- 
lution and the possible modifications discussed so far rely on two 
important assumptions which have recently been cast into doubt. 
First of all, the models considered in this work assume a simple re- 
lation between QSO activity and the mass of its hosting halo. How- 
ever, [w^jth^^^ogiJ J2005bll pointed out that the tight relation be- 
tween the BH mass and the velocity dispersion of the spheroid im- 



plies that it is the spheroid rather than the halo which determines the 
growth of the SBHs and the subsequent QSO activity. As a conse- 
quence, the observed correlation between the halo and BH masses 
should not be regarded as fundamental as it merely reflects the fact 
that massive haloes preferentially host bulges with large velocity 
dispersions. This would imply that QSO activity, closely related 
to the evolution of bulges, should be studied using the more so- 
phisticated models for galaxy formation and evolution. The second 
and more important issue is related to the results of the recent Mil- 
lennium Simulation fcpringel et alj2005l) The analyses performed 
bv lGaoetail 120051) and bv lHarker et alj <2005t) have shown that 
the spatial correlation properties and the formation epochs of the 
haloes depend on the local overdensity. This effect, which is partic- 
ularly evident for galaxy-sized haloes, contradicts one of the basic 
assumption of the EPS theory which we have used to construct the 
halo merger history and to describe their biasing function. It is not 
clear, however, how serious the implications are for galaxy/QSO 
formation models and for halo models of clustering. In case they 
are and in absence of a generalized EPS theory capabl e of account- 
ing for environm ental effects (see, however. I Abbas & Shethl2005t 
IShen et aljEoolh . then the only way out would be that of resort- 
ing to halo merger histories extracted from numerical experiments 
that implicitly account for environmental dependencies (see, e.g. 
iLemson & Kauffmanrll999T) . We will pursue this strategy in a sub- 
sequent work. 
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